Measurement of the Neutrino Asymmetry Parameter B in Neutron Decay 
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A new measurement of the neutrino asymmetry parameter B in neutron decay, the angular correla- 
tion between neutron spin and anti-neutrino momentum, is presented. The result, B = 0.9802(50), 
confirms earlier measurements but features considerably smaller corrections. It agrees with the 
Standard Model expectation and permits updated tests on "new physics" in neutron decay. 
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Assuming only Lorentz invariance, the decay probabil- 
ity for polarized neutrons is given by [1] 
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Pe, Pu, E, E, are momentum and energy of electron and 
anti- neutrino (in the following called neutrino), respec- 
tively, me is the electron mass, (s„) is the neutron spin 
polarization, and the fi^ denote solid angles. The param- 
eters a. A, B, and D are angular correlation coefficients: 
a is the correlation between the momenta of electron and 
neutrino. The parity violating parameters A and B for 
electron and neutrino asymmetry correlate the neutron 
spin with the momentum of electron and neutrino, re- 
spectively. The subject of this paper is a precise mea- 
surement of B. 

The factor K and the correlation coefficients are re- 
lated to the couplings of the theory. Within the Stan- 
dard Model of Particle Physics (SM), h vanishes and 
D is much smaller than the present experimental limit. 
K = Gl\Vud\'^F{E){g'^ + ?,g\), with Fermi-constant Gp, 
quark mixing matrix element Vudi phasespace factor 
F{E), and vector and axial- vector coupling constants 
gv and gA- AH correlation coefficients are functions of 
\ = gA/ gv ^ where gA and gv are assumed to be real, e.g.: 
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Small higher order corrections must be considered addi- 
tionally. The structure of the weak interaction {V — A 
in the SM) is not predicted by theory but has to be de- 
termined experimentally. Due to its sensitivity to the 
neutrino helicity the neutrino asymmetry parameter B 
is an important input parameter for this purpose: Pre- 
cise measurements oi B 0, 0] can be used to derive 
limits on hypothetical right-handed current [V+A) con- 
tributions to /3-decay [a H Q (mediated by new heavy 



charged bosons Wr)- These can be compared to other 
indirect measurements in muon [§] and nuclear /3-decay 
as well as to collider attempts to directly produce 
the Wr (cf. also [uj). Whe n interpreted in general 
left-right symmetric (LRS) models beyond manifest LRS 
theory, information from all these experiments (direct, 
muon, and /3-decay) are complementary and necessary 
to obtain limits [9|. 

A deviation from eq. ([2]) is a signal for new physics 
not described within the Standard Model and may be 
due to admixtures of right-handed currents, or due to 
anomalous (scalar, tensor) couplings pos sibly caused by 
exotic models like leptoquark exchange |12| . A recent 
review on this topic can be found in [l3| . Neutron decay, 
involving all particles of the first generation, is well suited 
to study the structure of the weak interaction [l^l since 
theoretical corrections are small and well calculable as 
they do not depend on nuclear structure [l^ 3, 17 1. 

In order to measure the neutrino asymmetry param- 
eter B, the electron spectrometer PERKEO II [18[ was 
installed at the cold neutron beam position PFIB [l^ 
at the High Flux Reactor of the Institut Laue-Langevin 
(ILL), Grenoble. A cold neutron beam with a thermal 
equivalent flux of 1.3x10^" n cm^^ s^^ was transversally 
spin polarized in a system of two supermirror polariz- 
ers in the new X-SM geometry A radiofrequency 
(rf) spinflipper [2l| allowed to reverse the spin direc- 
tion. Polarization P and spinflipper efficiency F, for- 
merly sources of large corrections and uncertainties, were 
determined to P = 0.997(1) and F = 1.000(1) leading to 
a 0.30(14)% correction on B. P and F were measured as 
a function of the neutron wavelength by a time-of-flight 
method. A second rf-flip per and two supermirror ana- 
lyzers in the geometry of [20| were used to determine F 
and to scan P. For the precise measurement of the abso- 
lute beam polarization, several opaque polarized '^He cells 
2^ , 23 1 were employed. Both, P and F were uniform over 
the full neutron beam cross section for all wavelengths. 
The flipping ratio, a measure for PF , was monitored reg- 
ularly and stayed constant during beamtime. 
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FIG. 1: The spectrometer PERKEO 11: Polarized neutrons 
pass the setup, the magnetic field divides the full solid an- 
gle into two hemispheres - in and against spin direction - 
and guides the decay products onto the detectors. The low 
energetic protons are accelerated onto a thin carbon foil on 
negative potential, where they generate secondary electrons 
that can be detected by the scintillators. 



Behind the polarizers, a neutron shutter was installed 
for background measurements in the electron detectors. 
It was made of enriched ^LiF, as were the neutron beam 
coUimation orifices and the beamstop at the end of the 
installation, since generation of 7-radiation and fast neu- 
trons in ^LiF is suppressed by 10"^ 



2JI . Additional shield- 



ing was employed to reduce the remaining background. 

A certain fraction of neutrons decayed within the de- 
cay volume centered in the spectrometer. Its main part 
is a pair of superconducting coils in split pair configu- 
ration, generating a magnetic field with a maximum of 
1.03 T perpendicular to the neutron beam (fig. [1]). The 
neutron spin aligns with the field that therefore separates 
the full solid angle in two hemispheres: One in and one 
against neutron spin direction. It guides the charged de- 
cay products onto the two detectors installed next to the 
beam, realizing a 2 x 27r detector where no solid angle cor- 
rections have to be applied. Systematic effects related to 
the spectrometer design are described below. A detailed 
description of a previous experiment can be found in [J] . 

The magnetic field B' decreases towards the detectors. 
This causes an increase of the parallel momentum compo- 
nent of the particles leading to reduced electron backscat- 
tering. This is further reduced by the "magnetic mirror 
effect" : Electrons scattered out of the detector may be 
reflected at the increasing B' and still detected in the 
correct hemisphere. Backscattering is recognized via its 
delayed signal in the second detector. The full energy 
of the decay electron is reconstructed. The fraction of 
events assigned to the wrong detector is smaller than 
0.2%, a neglectable systematic effect (< 10~*) if the re- 
gion of interest is chosen above a /3-energy of 240 keV 



Since the neutrino cannot be measured directly, elec- 
tron and proton were detected in coincidence to recon- 
struct the neutrino. Most sensitive to B is the case when 
electron and proton are emitted into the same hemisphere 
relative to the neutron spin - momentum conservation 
then restricts the neutrino to the opposite direction. The 
other case, where electron and proton are emitted into 
different hemispheres, is kinematically favored but less 
sensitive to B since the neutrino direction is less con- 
strained [2^ . As it depends strongly on detector calibra- 
tion, this case was only used for cross checks of the result 
obtained for the first case. 

Electrons (i?max — 782 keV) are detected by 
190 X 130 mm^ plastic scintillators with photomulti- 
plier readout. The protons having much lower energies 
(-^max — 780 eV) are accelerated onto a thin carbon foil 
(15—30 ^g cm^^, coated with MgO) on negative poten- 
tial (V^ = — 18 kV). Whereas the electrons pass the foil 
almost unperturbed, the heavy protons have enough ion- 
ization power to release one or several secondary electrons 
from the foil [23|. These are detected in the scintillator, 
where also the proton time-of-flight is registered. No pre- 
cise energy information on the proton is obtained with 
this method. Proton detection does not depend on the 
initial proton energy and the angles of incidence occur- 
ring in this setup as was experimentally verified. 

The measured signature is the experimental neutrino 
asymmetry defined by 
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where N^\E) is the number of coincident events with 
electron kinetic energy E. The first/second sign indi- 
cates whether the electron/proton was emitted in (-f ) or 
against (— ) neutron spin direction. Eq. ([3]) is related to 
the neutrino asymmetry parameter B by integrating eq. 
([1]) over the hemisphere [3, [2^ 
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The definition is separated into two regions by the energy 
dependent parameter r = P(E + 111^) / {Eynax — E) which is 
unity at £' = 236 keV. (3 = v/c. Eq. (|3|) is very sensitive 
to the coefficient B but also depends slightly on the cor- 
relations a and A whose experimental uncertainties have 
to be considered. 

Detector calibration was performed regularly and two- 
dimensional detector scans were carried out to correct for 
spatial detector characteristics. Due to the flat spectral 
shape of i?exp {E) detector calibration imposes only a tiny 
uncertainty of 0.02% on B. 

At low electron energies, there is background related to 
the high voltage (HV) applied on the carbon foils. Above 
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230-240 keV, however, the measured electron spectra, i.e. 
N^'^{E) and N~~{E), can be weh described by their 
theoretical expressions, where all fits have only one free 
parameter, a normalizing factor. An upper limit on re- 
maining background contributions in the fit region was 
determined from the fit residuals. At lower energies, a 
satisfactory description is impossible due to background, 
a non-linear energy calibration, and backscattered elec- 
trons assigned to the wrong detector [28| . 

All corrections due to the spectrometer design have 
been obtained from Monte Carlo simulations. The "edge 
effect" accounts for the loss of charged particles due to 
the finite length of the decay volume that was defined 
by thick aluminum baffles. "Grid Effect": Four layers of 
grids made of AlSi-wires (10 /zm and 25 /xm) were used to 
prevent the HV applied to the detector foils to reach into 
the decay volume. Different methods (finite elements, 
boundary elements) showed that the absolute electric po- 
tential in the decay volume is at least one order of mag- 
nitude below a value that would cause systematic effects 
at the present level of experimental precision. However, 
electrons and protons may be absorbed or scattered by 
the grids. This "grid effect" was obtained using the pro- 
gram PENELOPE [29,'| to simulate the electron trajecto- 
ries in the wires. Protons hitting the wires were assumed 
to be absorbed. The potential barrier for electrons to 
reach the scintillator can be neglected since all electrons 
with £' > 84 keV will certainly pass it regardless of their 
initial emission direction. 

Charged particles moving in an increasing magnetic 
field B' may be reflected as p\/B' is an adiabatic in- 
variant, where p± is the momentum component perpen- 
dicular to the field lines. This gives rise to the "mag- 
netic mirror effect" since a certain fraction of decay prod- 
ucts was emitted towards the field maximum due to 
the finite neutron beam width. An asymmetric setup, 
i.e. a displacement A between neutron beam and mag- 
netic field maximum, may cause an additional, possibly 
large effect on B. Therefore A was measured directly 
and was additionally determined from data in two inde- 
pendent ways to correct for the effect: It was obtained 
from a x^-minimization of fits to the difference spectrum 
D{E) = N—{E) - N++{E) that has the highest sensi- 
tivity on A at high electron energies E. A was also deter- 
mined from the relative difference of the electron asym- 
metry parameters A measured without ep-coincidence 
with the two detectors. This was possible due to the 
symmetric setup and since beam related background, i.e. 
background generated in the coUimation system that can- 
not be measured separately, is small (< 10"'*) and also 
symmetric in a region E > 350 keV [soj. The result- 
ing values for A do virtually not depend on B nor on 
other systematic effects. The associated error of 0.32% 
constitutes the largest systematic uncertainty of the mea- 
surement. 

F{E) in eq. ^ has already been corrected for Coulomb 



interactions Fc{E), proton recoil R{E), and outer radia- 
tive corrections Sii{E), and reads 

F{E) = F'iE) (1 + SniE)) (1 + R{E)) FdE), (5) 

where F'{E) is the uncorrected phasespace factor. The 
expressions for 6r{E) and RiE) were taken from [l^ and 
[3l|, respectively. The recoil and order-a corrections to 
B are of order 0.01% [11]. 

The proton coincidence window Wi was restricted to 
40 lis causing a small correction of —0.05(3)% to account 
for slower protons. Accidental coincidences were directly 
measured with a delayed coincidence technique in a de- 
layed window W2 from 42—82 ^s after the initial trig- 
ger. In order to avoid suppression of protons by acciden- 
tal coincidences or background, or suppression of acci- 
dental coincidences by preceding signals, up to 32 stops 
were detected in both detectors. Only events with ex- 
actly one stop in the respective window were considered 
in the analysis since multiple stops are mostly due to 
background. Events with a stop in Wi and a second 
("accidental") stop in W2 were included in the analysis. 

However, this "1 stop" condition causes an overesti- 
mation of accidental coincidences since the stop-signal 
combination "proton and accidental signal" in Wi is re- 
moved from the data set, whereas a similar combination 
does not occur in W2. The available information on all 
stops allowed to determine the necessary correction di- 
rectly from the data. 

At high electron energies E, the fit region is limited as 
the uncertainty related to the displacement A between 
neutron beam and magnetic field increases with E. At 
the low energy side, the fit region is limited by the ef- 
fects mentioned above: Background, non-linear detector 
response, and wrongly assigned backscatter events. The 



TABLE I: Neutrino asymmetry B: corrections and errors 



Detector 1 Detector 2 



Effect [%] 


Corr. 


Err. 


Corr. 


Err. 


Polarization 


+0.30 


0.10 


+0.30 


0.10 


Flip Efficiency 




0.10 




0.10 


Data Set: Statistics 




1.22 




0.36 


Proton Window 


-0.05 


0.03 


-0.05 


0.03 


1 Stop Condition 


-0.24 


0.06 


-0.13 


0.03 


Background 




0.10 




0.08 


Detector Calibration 




0.02 




0.02 


Spectrometer: Edge Effect 


-0.16 


0.05 


-0.16 


0.05 


Grid Effect 


-h0.03 


0.05 


+0.03 


0.05 


Mirror Effect 


+0.44 


0.05 


+0.44 


0.05 


Displacement A 


-0.10 


0.32 


+0.10 


0.32 


Correlations A, a 




0.07 




0.07 


Sum 


+0.22 


1.28 


+0.53 


0.52 
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FIG. 2: Fit of Bcxp to all detector 2 data. The solid curve 
indicates the fit region. The result does not depend on this 
region. However, the overall uncertainty increases if the fit is 
extended to higher energies due to the magnetic mirror effect. 



apply this value to analyze a manifest LRS model with 
zero mixing = 0), following the procedure described 
in [3] but with only the electron asymmetry parameter 
A from [m as further input parameter; A is a free pa- 
rameter as it may be different from the SM value. We 
obtain a lower Hmit mw^ > 290.7 GeV/c^ (90% CL). 
The rather small improvement of this limit despite the 
reduced uncertainty of B originates from the shift of B to 
a lower value. Due to the controversial neutron lifetime 
(cf. [lH) we renounce a more elaborated analysis. How- 
ever, with this controversy being settled, the improved 
neutrino asymmetry parameter together with other neu- 
tron decay data will permit to derive new limits on gen- 
eral LRS models and on scalar and tensor couplings. 



region was chosen from 250—455 keV. The final asymme- 
try parameter B is independent of this choice as the fit 
results agree within ib0.3gst,a,t, for intervals between 235 
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and 620 keV (fig. 4.48 in (281). 

Fig. [5] shows the fit of i?cxp, eq. to all data of de- 
tector 2 (proton efficiency about 17%). ft yields the neu- 
trino asymmetry parameter B2 = 0.9798(36)stat(36)syst. 
The result of detector 1, Bi = 0.9845(120)stat(36)syst, is 
limited by statistics due to a smaller proton efRciency. 
This was caused by an inferior foil coating and higher 
HV background that could not be further suppressed. A 
detailed compilation of all relevant corrections and errors 
is given in table [D 

In this situation, with two detectors of very different 
statistical significance, we use the statistical average as 
the final neutrino asymmetry parameter result 



B ^ 0.9802(50) = 0.9802(34)stat(36) 



syst- 



(6) 



The experiment is limited by statistics and the uncer- 
tainty due to the displacement A between neutron beam 
and magnetic field. With two detectors of equal perfor- 
mance, both errors would be significantly smaller as the 
influence of A would cancel by calculating the arithmetic 
mean of Bi and B2 ■ 

The second case, electron and proton in opposite detec- 
tors, yields results with much larger uncertainties, 1.9% 
and 3.0% for detector 1 and 2, respectively, dominated 
by detector calibration. They statistically agree with ([6]). 

Our result ([6]) has a precision similar to the most pre- 
cise measurement so far Q and agrees very well with all 
results published earlier. It is distinguished, however, as 
it features several times smaller corrections than com- 
peting experiments (0.5%; 1% if absolute numbers are 
considered, cf. table |T|). Our result is consistent with 
the Standard Model expectation, Bsm = 0.9876(2), cal- 
culated with the current world average for A from [lll| 
and eq. 

By including our result (O, the error of the world aver- 
age for B reduces by 25%, yielding B = 0.9807(30). We 



[1 
[2: 

[s: 
[4: 

[5] 
[6] 
[7 

[s: 
[9: 
[10: 

[11 

[12 
[13 

[14 

[15 
[16 
[17 
[18 
[19 
[20 
[21 
[22' 
[23^ 
[24' 
[25^ 
[26' 
[27 

[28; 

[29 
[30 

[31] 



* Electronic address: |marc.schumann@gmx.net|_ 



Electronic address: abele@physi.uni-heidelber g.de| 
J. D. Jackson et al., Phys. Rev. 106, 517 (1957) 
I. A. Kuznetsov et al., Phys. Rev. Lett 75, 794 (1995) 
A. Serebrov et al., JETP 86, 1074 (1998) 
M. Kreuz et al., Phys. Lett. B 619, 263 (2005) 

D. Dubbers, Prog. Part. Nucl. Phys. 26, 173 (1991) 

A. Serebrov, N. Romanenko, JETP Lett. 55, 503 (1992) 

H. Abele, Nucl. Instr. Meth. A 440, 499 (2000) 

J. R. Musser et al., Phys. Rev. Lett. 94, 101805 (2005) 

E. Thomas et al., Nucl. Phys. A 694, 559 (2001) 

V. M. Abazov et al., Phys. Rev. D 69, llllOl(R) (2004) 

W.-M. Yao et al. (PDG), J. Phys. G 33, 1 (2006) 

P. Herczeg, Prog. Part. Nucl. Phys 46, 413 (2001) 

N. Severijns, M. Beck, O. Naviliat-Cuncic, Rev. Mod. 

Phys. 78, 991 (2006) 

H. Abele, Part. Prog. Nucl. Phys., 

doi:10.1016/j.ppnp.2007.05.002 (2007) 

A. Sirhn, Phys. Rev. 164, 1767 (1967) 

F. Gliick, Phys. Lett. B 436, 25 (1998) 

S. Ando et al., Phys. Lett. B 595, 250 (2004) 

H. Abele et al., Phys. Lett. B 407, 212 (1997) 

H. Abele et al., Nucl. Instr. Meth. A 562, 407 (2006) 

M. Kreuz et al., Nucl. Instr. Meth. A 547, 583 (2005) 

A. N. Bazhenov et al., Nucl. Instr. Meth. 332, 534 (1993) 

O. Zimmer et al., Phys. Lett B 455, 62 (1999) 

O. Zimmer, Phys. Lett. B 461, 307 (1999) 

M. A. Lone et al., Nucl. Instr. Meth. 174, 521 (1980) 

M. Schumann, H. Abele, arXiv:0708.3150 (2007) 

F. Gliick, I. J06, J. Last, Nucl. Phys. A 593, 125 (1995) 

D. E. Kraus, F. A. White, IEEE Trans. Nucl. Sci. NS-13, 

765 (1966) 

M. Schumann, PhD thesis. University of Heidelberg 

www . ub . uni-heidelberg . de/archiv/7357 (2007) 

J. Sempau et al., Nucl. Instr. Meth. B 207, 107 (2003) 

D. Mund, PhD thesis. University of Heidelberg, 

www . ub . uni-heidelberg . de/archiv/ 6576 (2006) 

D. H. Wilkinson, Nucl. Phys. A 377, 474 (1982) 



